Probing level renormalization by sequential transport through double quantum dots 
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We study electron transport through double quantum dots in series. The tunnel coupling of the 
discrete dot levels to external leads causes a shift of their energy. This energy renormalization affects 
the transport characteristics even in the limit of weak dot-lead coupling, when sequential transport 
dominates. We propose an experimental setup which reveals the renormalization effects in either 
the current-voltage characteristics or in the stability diagram. 
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I. INTRODUCTION 

Serial double quantum dots are ideal systems to inves- 
tigate various quantum mechanical effects such as molec- 
ular binding 1,2 or coherent dynamics 3 between the con- 
stituent dots. Furthermore, they are considered as an 
implementation of a charge 4 or spin qubit. 5 Elaborate 
experimental techniques were developed to control and 
characterize double-dot structures, 6-8 and many informa- 
tion about the system can be deduced from the electric 
conductance through the device. 9 Recent experiments in- 
clude the measurements of quantum mechanical level re- 
pulsion due to interdot coupling 10 as well as due to ex- 
ternal magnetic fields, 11 the detection of molecular states 
in a double dot dimer, 12 and the observation of coherent 
time evolution of the dot states. 3 

Transport through serial double dots, as depicted in 
Fig. 1, inherently visualizes the basic quantum mechani- 
cal concept of coherent superposition of charge states. 13 
The states that are coupled to the left and right lead, 
the localized states in the left and right dot, respectively, 
are no energy eigenstates of the double dot. This leads 
to oscillations of the electron in the double dot as it was 
shown in recent experiments. 3 ' 4 To account for this in- 
ternal dynamics, descriptions using classical rates only, 
are insufficient, which is why approaches including non- 
diagonal density matrix elements for the double dot have 
been developed. 14-17,23 . 

In this paper, we propose to use a serial double quan- 
tum dot to probe another consequence of quantum me- 
chanics: the energy level renormalization of the quantum 
dot levels due to tunnel coupling to a reservoir. This idea 
is based on two properties of a serial dot system. First, 
the left and right dot levels are tunnel-coupled to different 
reservoirs. Since the level renormalization is a function 
of the level energy, the lead chemical potential, and the 
tunnel coupling, the energy shift of the left and right dot 
levels is, in general, different. Second, the conductance 
through the double dot is very sensitive to the difference 
of the energy levels. It shows a resonant behavior with 
the width given by the tunnel couplings, 18 which can be 
much smaller than the temperature. This sharpness of 
the resonance makes the conductance a valuable experi- 
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FIG. 1: Schematic energy profile for a double dot coupled in 
series to two reservoirs. Each reservoir is coupled to the dot 
of the corresponding side by the coupling strength F r . The 
inter-dot coupling is determined by A. The energies of the 
two dot states are characterized by the mean energy E and 
their relative distance e. 



mental tool, for example to measure the shell structure 
of quantum dots. 19 

It is well known 20,21 that tunnel coupling to reservoirs 
renormalizes the energy levels. In single-dot geometries 
such an energy renormalization is only accessible in trans- 
port of higher order in the tunnel coupling strength. As 
we will show below, this is different for the serial dou- 
ble dot geometry, for which renormalization effects are 
visible in the conductance already in the limit of weak 
dot-lead coupling, described by transport to first order 
(sequential tunneling) in the tunnel- coupling strength 

r = r L + r R , 

The paper is organized as follows: In Sec. II we present 
the model Hamiltonian for the double dot. 15-17,23,23 ' 24 
In Sec. Ill we derive the stationary density matrix and 
the dc— current for arbitrary bias voltages. In Sec. IV, 
we discuss how renormalization effects appear in the 
dc— current through the double dot. We draw our con- 
clusions in Sec. V, followed by some technical notes in 
the appendices. In App. A we make explanatory notes of 
the diagrammatic technique, and in App. B we give an 
illustrative reformulation of the master equation in terms 
of a pseudo spin. 
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II. MODEL 

We consider a double quantum dot, contacted in series, 
which is described by the Hamiltonian 23 

H = Hr + Hd + H ? ■ 

r=L,R 

The first part of the Hamiltonian describes the electric 
contacts on the left (L) and right (R) side. These con- 
tacts are modeled by large reservoirs of noninteracting 

electrons H r = Tr k ,a e rk c tka c rka- Hcrc c rk<r,4ko denote 
the annihilation and creation operators for electrons in 
the reservoir r £ {L,R} with spin a. The reservoirs are 
assumed to be in equilibrium, so that they can be char- 
acterized by the Fermi distribution /l/r(w). An applied 
bias voltage V is modeled by different chemical potentials 
in the left and right contact /l/r(w) = f(w ± eV/2). 

The second part of the Hamiltonian, Hy>, describes 
two dots, containing one electronic level each, which are 
coupled by the Coulomb interaction: 

H d = Eririr + Un^n^ + U 1 {n^ni^i+n^n^) . (2) 

r=L,R 

Here, n ia = df a d ia and = df a d irT are the occupa- 
tion number operators for dot i £ {L,R} with spin tr, 
where di,y,df a being the annihilation and creation oper- 
ators of an electron on dot i with spin a. Each dot con- 
sists of a single electronic level at the energy E^ /r mea- 
sured relative to the equilibrium chemical potential of the 
leads. We parameterize the levels by their average energy 
E = (Eh + Er)/2 and their difference e = El — £r, so 
that £l/r = E±e/2. Double occupation of one individ- 
ual dot is associated with the intradot charging energy 
V. Simultaneous occupation of the both dots with one 
electron each costs the interdot charging energy U with 
U' 3> U. States with three or more electrons in the dou- 
ble dot are not considered in the following. The remain- 
ing eigenstates of Hy>, then, are: both dots empty |0), one 
electron with spin a in the left |Lct) or right dot \Ra), and 
one electron in each dot |LctR(t'). We assume that the 
intra-dot charging energy always exceeds the lead Fermi 
energies. Therefore the states with two electrons in the 
same dot |LerL<r) and IRctRct) will have a vanishing occu- 
pation probability. However, these states will appear as 
intermediate (virtual) states in our calculation, providing 
a natural high-energy cut-off. 

The third part Ht = Ha + Hy of the Hamiltonian 
Eq. (1) describes both, tunneling between the two dots, 
Ha, as well as tunneling between dots and leads, Hy, 

ffA = -yE(<^ + 4.Aj (3) 

H r = Y1 *Lfc c Lfc<A CT + tnkORkadaa + h - c - ( 4 ) 

ka 

Due to the serial geometry, an electron from the right 
(left) reservoir can only tunnel to the right (left) dot. 



The tunnel coupling of reservoir r to the corresponding 
dot is characterized by the coupling strength T r (u>) = 
^ n ^2k \trk\ 2 ${£rk — w). We consider only spin conserv- 
ing tunneling processes, and assume flat bands in the 
reservoirs, which yields energy independent couplings T r . 
Furthermore, we choose the interdot tunnel coupling am- 
plitude A as a positive, real parameter, which can be 
always achieved by a proper gauge transformation. 

III. KINETIC EQUATION 

In the following section we calculate the stationary re- 
duced density matrix p s t for the double dot system and 
the dc— current through the system. The reduced density 
matrix of the double-dot is obtained from the density ma- 
trix of the whole system by integrating out the reservoir 
degrees of freedom. The Liouville equation for the re- 
duced density matrix then has the following structure: 

= ih-^Pst = [Hjj, pst] + [Ha, Pst] + Sp st • (5) 

The first two parts represent the internal dynamics on 
the double dot, which depends on the level separation e 
and the interdot coupling A. The third part of Eq. (5) 
accounts for the tunnel coupling between double dot and 
external reservoirs. The fourth order tensor S contains 
imaginary and real parts, associated with particle trans- 
fer processes and with tunnel induced energy renormal- 
ization of the dot levels, respectively. The latter has been 
neglected in previous works. 16 ' 17,23 We calculate E using 
a real-time diagrammatic approach 21,25 as explained in 
App. A. Also alternative methods are available such as 
Bloch-Rcdficld theory. 22,26 

In the following wc concentrate on the limit of weak 
tunnel coupling between double dot and leads. There- 
fore, we calculate S to lowest order in the tunnel- coupling 
strength T = Tl + Tr, which defines the so-called 
sequential-tunneling approximation. This approximation 
implies that all tunneling events are independent from 
each other, which is fulfilled for /cbT ^ T. Since the 
correlations generated in the bath during a tunnel pro- 
cess decay on the time scale h/k^T, 27 (this follows from 
the dependence of the tunneling line in Fig. 4-6 on its 
extension in time), while the average time between con- 
secutive tunneling events is given by the inverse of the 
coupling strength h/T, higher order, coherent tunneling 
events are suppressed by the condition k^T T and 
may be neglected. 

The energy eigenstates of the double dot subsystem 
Hy> + Ha are the bonding and anti-bonding states with 
energies E^u= E=p A a b where A a b = VA 2 + e 2 denotes 
their energy splitting. This identifies A a b as frequency of 
the charge oscillations, 3,4 and A as minimum distance be- 
tween the bonding and anti-bonding eigenstates as func- 
tion of the left and right energy level. 10 

If the splitting exceeds the intrinsic broadening of the 
levels, A a b ^> T, then the internal oscillations are fast, 
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and transport through the double-dot system takes place 
through two separate incoherent levels. In this case, off- 
diagonal matrix elements of the stationary density ma- 
trix vanish, which can be seen from the expansion of the 
Liouville equation = ih-^P^ = A ah P£ + 0(T), where 
P* 1 denotes the matrix element P*^ 1 = (xilPst IX2) of 
the reduced density. 

The more interesting transport regime is in the oppo- 
site limit, A a b < r, where the external coupling strongly 
modifies the internal dynamics, which is captured by the 
off-diagonal elements of the reduced density matrix. 17,23 
Combined with the validity condition for sequential tun- 
neling, i.e. r <C k^T this implies A a b -C k^T, i.e., 
internal oscillations are slow as compared to the time 
scale for the correlations during a tunneling event. As a 
consequence, the localized states |L<r) and |Rcr) can be 
used as eigenstates of the double dot in the calculation 
of E, which facilitates the interpretation of the dynam- 
ics. Technically, the condition A a b < T <C k-gT means 
that for a consistent theory, we do not only expand £ 
in Eq. (5) to first order in T, but also have to expand it 
to zeroth order in A a b- This is accomplished by replac- 
ing the energies Pl/r arising in the calculation of £ by 
the mean level energy E = (Pl + Pr)/2. (Therefore our 
formulas only contain the Fermi functions at the average 
single particle level / r (P), while energies of the order of 
the interdot tunneling or the level separation are smeared 
by temperature). It is worth to point out, that by using 
the localized states as basis of the £ (i.e. by calculating 
the transition rates in the localized basis) one is auto- 
matically limited to the regime A, e < fc B P. 14,16,17,23 

The technical details of how to calculate £ are de- 
scribed in App. A. The master Eq. (5) then must be 
solved under the constraint of probability normalization 
Tr[p s t] = 1. The stationary current / is given by the time 
derivative of the expectation value of the total number 
of electrons in either the left or the right lead. For the 
lowest-order expansion used in the present context, the 
current can alternatively be written in the form 23 



I = -e^(lH A ,n L ])=-^AIm(J2P&) , (6) 



where Im denotes the imaginary part. In App. B we give 
an analytical solution for the current as function of bias 
voltage and gate voltages. In the following we discuss 
our results. 

Instead of working with an off-diagonal density matrix 
one may switch to a pseudo spin representation of the 
problem. Then this double-dot transport problem shows 
similarities to the system of a quantum dot connected to 
ferromagnetic leads. 28,29 This will be discussed in App. B. 



IV. DISCUSSION 

The stationary current takes the form 

/ = -A 2 - (7) 

h P 2 + £ 2 cn [l) 

The numerical factors A and B (the explicit form is given 
in App. B) depend only on the tunnel coupling constants 
Tl, Tr, and A as well as on the Fermi distribution func- 
tions /l/r(P) and /l/r(P + U) of the left and right lead, 
but not on the level energy difference e. The current in 
Eq. (7) shows the well-known 18,23 Lorentzian dependence 
on the energy separation e ren between left and right dot 
level. However the energy separation is affected by the 
renormalization of the bare localized levels: 

e rcn = e + AE L - AP R . (8) 

This is the central statement of the present paper. The 
energy shift AE r of the energy level in dot r, caused by 
the external tunnel coupling, is given by 

AE r = <j) r (E)~2<j) r (E + U) + (j) r {E + U') (9) 

with 

^H^Re*^^). (10) 

Here, Re denotes the real part, 4" is the digamma func- 
tion, and ^l/r = ±eV/2 the leads' chemical potentials. 
We want to emphasize that this energy renormalization 
is not due to a capacitative but rather due to the tun- 
nel coupling to the reservoirs. Furthermore, it vanishes 
in the noninteracting case U = U' = 0. The intradot 
charging energy U' (which we usually treat as infinite to 
avoid double occupation of one dot) serves as a natural 
cut off for the energy renormalization in Eq. (9). This is 
the reason why we allowed the intermediate states \5 m 
App. A to occupy these states. 

The energy shift of the localized levels is proportional 
to the tunnel coupling strength and depends on the dot 
level positions relative to the Fermi energy. The renor- 
malized level separation as function of the bias voltage is 
plotted in Fig. 2a. The renormalized level separation e rcn 
reaches a (local) extremum each time, when the Fermi en- 
ergy of a lead becomes resonant with the energy needed 
for single (/j. r = P) or double occupation (/j, r = E + U). 

Fig. 2b shows the current as function of the transport 
voltage taking the level shift into account (solid line). By 
neglecting the level shifts (dashed line in Fig. 2b), the 
current shows a typical Coulomb staircase. The steps 
occur when a lead chemical potential aligns with an elec- 
tronic level in the double dot. Since the bare energy level 
separation e as well as the interdot tunneling A shall be 
of the order of or smaller T > {A, e} and we consider 
r < /cbT, the different single particle states are not re- 
solved as individual steps in the I — V staircase. The 
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FIG. 2: Upper Part: Renormalized level spacing e rG n (solid 
line) between the electronic levels in the left and right dot as 
function of the transport voltage V. e rGn is extremal, when 
the chemical potential of a lead aligns with the energy needed 
for either single (E) or double occupation (E + U). Lower 
part: Current- voltage characteristics for bare (dashed line) 
and renormalized level spacing (solid line). Rcnormalization 
of energy levels leads to an asymmetric current- voltage char- 
acteristic. The current increases (decreases) whenever the 
level spacing is reduced (increased) with respect to the bare 
value. Plot parameters are: e = A = Tr = Tl = T/2, 
E = 10/cbT, U = 20fc B T, and U' = 100k B T. 



tunnel induced renormalization leads to additional fea- 
tures on the staircase (solid line in Fig. 2b). Whenever 
the magnitude of the renormalized level spacing grows 
(drops) the current decreases (increases). This leads to 
a suppression or an enhancement of the current around 
the steps of the I — V characteristic, leading to regions 
of negative differential conductance. The width of these 
feature is of the order of the charging energy and can 
exceed temperature and coupling strength significantly. 

Neglecting renormalization effects and assuming sym- 
metric coupling to the reservoirs (r^ = T^), the current 
through the double dot is an odd function of the trans- 
port voltage (see dashed line in Fig. 2b). This is no longer 
the case when renormalization is taken into account (see 
solid line in Fig. 2b). The reason for this asymmetry is 
that even though the change of asymmetry, AEl — AEr,, 
caused by level rcnormalization is antisymmetric with re- 



spect to the bias voltage, this in not true for the total 
asymmetry e ren = e + A£l — ASr due to the non- 
vanishing bare splitting e (see Fig. 2a). A comparable 
asymmetry in transport through two coupled dots was 
recently observed by fshibashi et al. 30 and theoretically 
described by Fransson et al. 20 However, a negative dif- 
ferential conductance feature can not be uniquely linked 
to such renormalization effects. Due to interface capaci- 
ties the level positions in the left and right dot are always 
affected by the transport voltage in real experiments. 18 ' 31 

To exclude the effect of interface capacities, we propose 
a different experiment: measuring the current I{E\ J1 Er) 
at a constant transport voltage as function of the left 
and right gate voltages on the dots. The resulting sta- 
bility diagram is plotted in Fig. 3a. Elastic sequential 
tunneling from the left to the right dot is possible if 
El ~ E-r. Furthermore electron transport from the left 
to the right reservoir takes only place if the dot level for 
single (E) or double occupation (E + U) is located in the 
bias voltage window. Therefore the current resonance 
forms two stripes in the regions —eV/2 < E < eV/2 and 
-U - eV/2 < E < -U + eV/2. Away from the current 
stripes the occupation number of the left and right dot 
(JVl, JVr) is fixed, and no current can cross the structure. 
For a detailed discussion on stability diagrams for trans- 
port through double dots we refer to the review of van 
dcr Wicl. et al. 6 

In the absence of renormalization effects, the cur- 
rent stripes would exactly coincide with the condition 
£l = En. By plotting the current as function of the 
mean level position E = (E^ + E-r)/2 and the relative 
energy difference e = £l — -Eft, one would therefore ex- 
pect a straight horizontal line. Instead, the maximum of 
the current follows the renormalization shift, where the 
condition e Icn = is fulfilled, sec Fig. 3b. The shift of 
the resonance is of order T as shown in Eq. (9) and can 
be small on the scale of bias voltage or temperature. The 
width of the current maxima in the stability diagram in 
Fig. 3 is not determined by temperature but rather by 
the dominant coupling strength 18 max(r, A). Therefore 
the resonance width is sharp enough to be able to mea- 
sure the renormalization of energy levels if T > A as used 
in Fig. 3. 

In the nonlinear transport regime AE r depends on E 
and therefore the current stripes in Fig. 3 are bent and 
tilted against each other. This dependence can be used 
as a stringent experimental prove of the renormalization 
of energy levels. Due to internal cross capacities, always 
appearing in real experiments, the gate voltage of one dot 
is a linear function of the gate voltage of the other dot. 
Therefore the stability diagram /(Vl,Vr) as plotted in 
Fig. 3a would experience a linear shear transformation. 
However straight (parallel) lines stay straight (parallel). 
Thus, cross capacities can not mimic the bending due to 
renormalization effects. 

In real experiments in addition to the resonant current 
stripes explained here, further features can arise due to 
inelastic processes, cotunneling, or due to excited levels 
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FIG. 3: Upper panel: Stability diagram I(El,Er) of the 
current through the double dot in the nonlinear transport 
regime. Well inside the areas, separated by the black line, the 
occupation of the individual dots is fixed to the written values 
(Nl,N-r). Elastic sequential current can cross the structure 
for E L « E-r and either -eV/2 < E < eV/2 or -eV/2 < 
E + U < eV/2 resulting in two current stripes. Lower panel: 
Gray scale plot of the current as function of the average single 
particle energy E and the bare level separation s = El. — E-r,. 
The different renormalization of left and right level shifts the 
current maxima by AEr, — AEl (dashed black-white line) 
where e rcn = . This leads to a tilting of the current stripes 
relative to each other. Relevant plot parameters are feT = 

2r, r L = r R , = a = 172, v = wk B T, u' = ioofc B r, and 

interdot charging energy U = 20kBT. 



within the bias voltage window. 6 These effects mainly 
lead to features within the triangles below the current 
strips in Fig. 3a and are expected not to interfere with 
our presented results. 

Finally we compare our result obtained for the sta- 
tionary current Eq. (B7) with previous theoretical works. 
For this we set the Fermi functions to fh(E) = 1 and 



I = 



f R (E + U) = 0. This simplifies 



r R A 2 



A 2 



In 
2r, 



+ 4(e rcn )2 + r| 



Neglecting renormalization effects (setting e rcn = e), 
this equation reproduces Eq. (4.19) in the paper by 
Gurvitz, 17 . Choosing the voltages such that the dot 
structure can also be doubly occupied, i.e. fh(E + U) = 
fh{E) = 1 and = f R (E) = f R (E + U) one obtains 
Eq. (4.18) of Ref. 17. 

Several publications assume, that if the lead Fermi en- 
ergies are far away from the electronic states of the dots, 
then the principal value integrals (Eq. B5), leading to 
the renormalization, can be neglected. However the en- 
ergy shifts are relevant on an energy scale given by the 
charging energy U, as shown in Fig. 2a. Therefore the 
assumption, that one can neglect renormalization effects 
and still exclude states with more than one electron oc- 
cupying the double dot is not justified. 



V. CONCLUSIONS 

If a quantum dot is connected to a reservoirs, the tun- 
nel coupling causes an energy renormalization of the elec- 
tronic states. We derived the conductance of a double 
dot connected in series to external reservoirs for general 
bias voltages and temperatures, taking into account these 
energy renormalizations. We have shown, that the con- 
ductance of such a double dot structure is affected by the 
energy level shifts already in a lowest order expansion in 
the tunnel coupling strength, due to its high sensitivity 
on the relative detuning of energy levels. Therefore we 
propose to use a double-dot system as detector for these 
energy renormalization effects. 

We present experimental consequences of the renor- 
malization in the current-voltage characteristics and in 
the stability diagram for the double dot in the nonlinear 
transport regime. In the current-voltage characteristics 
we find prominent negative differential conductances in 
voltage windows of the order of the charging energy. 

In the stability diagram of the double dot, we found 
that the current stripes arising as function of the gate 
voltages for left and right dot are tilted against each other 
and do not lie on a straight line, as it is the case when 
energy renormalization is neglected. We showed that the 
tilting of the current stripes is resolvable even in the se- 
quential tunneling regime (i.e. for F > k R T) as long as 
the interdot tunneling, A is of the same order or smaller 
than the external coupling T > A. 
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APPENDIX A: DIAGRAMMATIC RULES 
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In the following we show, how we calculate the tensor 
S X4X2' where Xi e {|0) , \La) , |Rer) , |LcrRcr')} arc the lo- 
calized eigenstates of H^, including the spin degree of 
freedom. 

We apply a diagrammatic technique, where each class 
of tunneling processes can be represented by a diagram. 
Its general derivation can be found in Ref. 21. Recently 
this technique was applied to Anderson-like Hamiltonians 
to investigate spin-valve effects, 28 or signatures of the 
excitation spectrum in the Coulomb blockade. 32 

Within this approach, the tensor is represented 

as block diagram, which is a part of the Keldysh time con- 
tour as shown in Fig. 4. The upper and lower line of the 
Keldysh time contour ij< represent the propagation of the 
double dot system forward and backward in time. They 
connect the matrix element characterized by the labels 
on the left side with the matrix element characterized by 
the labels on the right side. In the sequential tunneling 
approximation all transitions are allowed where a single 
electron first leaves and then reenters the double dot or 
vice versa. The two tunnel Hamiltonians are represented 
by vertices on the propagators. These vertices are con- 
nected by the contraction of the lead Fermi operators 
(indicated by a dashed line). Each line is characterized 
by its energy w, the spin a of the transfcred electron, as 
well as the corresponding reservoir label r 6 {L,R}. A 
vertex with an outgoing (incoming) tunneling line rep- 
resents an electron leaving (entering) the double dot on 
the specified side r. All possible transitions in lowest or- 
der in the external coupling T belong to one of the eight 
diagrams depicted in Fig. 5. 



VX3X1 

X4X2 



Xi' 



•X3- 



X2» -''' < X4 - ' 

^ ^real 



FIG. 4: Sketch of the structure of a diagram. The upper 
(lower) horizontal line denotes the forward (backward) prop- 
agator of the double dot system. The Keldysh time contour 
is labeled by tn, while the real time runs from left to right. 



FIG. 5: All topologically different diagrams contributing to 
the tensor ££4x2 ca l cu l a ted in first order in the external cou- 
pling r. Labeling of the eigenstates at the four corners and 
of the tunneling line like in first diagram, xs labels an inter- 
mediate charge state of the double dot. 



2^X4X2 is given by the sum of all diagrams with the 
corresponding eigenstates at the four corners, see Fig. 5. 
The number of relevant diagrams is limited by spin and 
particle number conservation as well as to the serial sys- 
tem geometry. The rules to evaluate these diagrams in 
lowest order are: 

1. Draw the upper and lower time contour. Add two 
tunnel vertices in any topological different way. 
The relevant criteria are the upper and lower con- 
tour, and the time ordering of the vertices on the 
real axes, not only on the Kelysh time contour. As- 
sign to each free segment of the contour a state of 
the double dot and the corresponding energy. For 
'bubble' diagrams like in the lower row of Fig. 5, 
an intermediate state %5 participates. 

2. The two vertices are connected by a tunnel line. 
Each tunnel line is labeled with the energy of the 
tunneling electron lu, its reservoir label r and its 
spin a. Spin and reservoir label of the tunneling 
electron are uniquely determined by the eigenstates 
involved in the tunneling processes. 

3. Assign to each diagram the resolvent 1 / {AE + i0 + ) 
where AE is the difference between energies be- 
longing to left going lines and energies belonging 
to right going lines (the tunneling line as well as 
the propagators). 

4. The tunneling line connecting two vertices and la- 
beled by the reservoir index r gives rise to the factor 

Here, the Fermi function ff(u>) = f r (u)) = 1/(1 + 
exp[(cj — ix r )/k^T\) corresponds to a tunneling line 
that is backward directed in the Keldysh time or- 
dering (compare Fig. 4), and /^(w) = 1 — f r (u>) 
corresponds to a tunneling line forward directed in 
the Keldysh time ordering. 

5. Each diagram gets a prefactor (— 1)° where v is 
the number of vertices on the backward propagator. 
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(This leads to an (—1) for the diagrams in the upper 
row of Fig. 5.) 

6. Sum over possible internal eigenstates %5 and inte- 
grate over the energy to of the tunneling electron. 
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FIG. 6: Relevant diagrams contributing to two specific entries 
of E, in a lowest order expansion in F Every diagram corre- 
sponding to a specific entry is labeled by the same eigenstates 
at its four corners. 

In the parameter regime we are interested in, the fol- 
lowing relations hold: kT > V > e, A. Therefore the 
energy difference between the single particle states is not 
resolved by the Fermi functions in the reservoir, so that 
we have to approximate the eigenenergies of {|0), |L<r), 
|Rct), |LcrRo-')} by {0, P L w E R w E, 2E + U}. While we 
exclude a double occupation of a single dot for the initial 
or final states by setting f r (E + U') = we allow the 
intermediate state X5 to be in such a state. These states 
have the eigenencrgy 2E + U' . 

In Fig. 6, we show as examples the diagrammatic ex- 
pansion of the tensor elements E^'J] and E^'^.. E^'q 
is purely imaginary and its magnitude has the meaning of 
a transition rate for a tunneling-in process starting from 
the empty double dot and resulting in a single electron 
with spin a sitting in the left dot. In contrast, E^'^ 
also has a real part which renormalizes the energy lev- 
els. Calculated in lowest order in T, each element of the 
tensor E can be expressed by terms of the form: 



7?H 



(A2) 



m(E-Lu) + iO+ ' 
where n and m are either (— ) or (+). In this notation, 



the algebraic expression for Ep^' 1 ^. is: 



j R<tR<t 



= x<r> + \E) + x$->-\E) 



(A3) 



+g a (x(+'+\E + U)+ Xp>-\E + [/)) 

+ (g a - 1) (xt+>+\E + U') + X^-\E + U'j) 

where, within this appendix, we allow for an arbitrary 
spin degeneracy g a . Since f r (E + U') = the imagi- 
nary part of the last row vanishes, however this is not 
the case for the real part, which causes the level renor- 
malization. The real part of the diagrams is determined 
by the principal values of the integrals in Eq. (A2) and 
can be expressed as a sum over digamma functions, see 
Eq. (9). 

Since the Hamiltonian given in Eq. (1) is indepen- 
dent of the orientation of the spin, each spin-realization 
of a charge state is equally probable. We can there- 
fore define P = (0|p st |0), P r r , = (ia\ p st \i'a), and 
P2 = J2a a> (LctRct'I p s t |LerRer'). Furthermore the sta- 
tionary density matrix is diagonal in spin and particle 
number. Thus the reduced density matrix p describing 
the double dot is given by the 4x4 matrix 



Pst 



/p 










pKpK 





pL 
R 







Pi 



\ 



(A4) 



The diagonal elements of the density matrix are the 
probabilities to find the double dot empty (Pq), the left 
(Pl ) or right dot (Pr) singly occupied, or the two dots 
simultaneously occupied by one electron (P2). Super- 
positions of the two single occupied states are possible 



pL 
R 



(p?r- 



One can define an effective tensor for E, that only de- 
pends on the orbital part of the matrix elements (denoted 
in the following formula by xi, X2, X3i Xa) and no longer 
on the spin variables. The new tensor elements are de- 
fined by: 



VX3X1 

X4X2 



V s x f 1 

XiX2 



(A5) 



Here i labels any possible spin-realization for the initial 
states, Xi,X2, and / for the final states X3iX4- (Due to 
spin degeneracies the two particle states are four fold de- 
generate, and the left and right states are each two-fold 
degenerate.) The tunnel tensor E^xi ^ s independent of 
the spin-realization i. The spin degeneracy appears only 
as a prefactor, but does not change the functional form of 
the elements. For example, E^'q = ^ CT eJ^'q describing 
the transition from Pq to Pl is twice as big for spin- 
degenerate electrons as for spin-less fermions. On the 
other hand = E$£} + E^j = E^ft + E^Jft 

describing the loss term of Pl is the same for spin- 

L I LT 

degenerate or spin-less fermions since E L j' L j = = 
v lt,U 
14,14 ■ 
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This treatment of the spin allows a general solution 
of the problem including both, the case of spin polar- 
ized electrons and the case of spin degenerate electrons. 
For the interested reader, we specify the degeneracy of 
fermions in the further Appendix by the variable g a : 
g a = 2 for electrons, g a = 1 for spin-less fermions. 



APPENDIX B: REWRITING KINETIC 
EQUATION AS BLOCH LIKE EQUATION 



by a Bloch-likc equation: 



dt 



(It 



-(Bxl) (B4) 



E A ^ [g°fr(E)P - 



1 



-{g a f r {E+U)-f-(E))P l 
f-(E + U)P 2 ] 

\Y,T^{E)+9.fr{E + U))l 



Instead of working with off-diagonal density matrix el- 
ements, we can switch to a pseudo spin representation. 
As any two level system, the 2x2 hermitian submatrix 
of the singly occupied states in Eq. (A4) can be treated 
as SU(2) representation of a pseudo spin Bloch vector 

1 = ( p k + P l\ iP k ~ iP hi p t ~ p r) T / 2 - For a complete 
set of variables, we further introduce Pi = P^ + P^ as 
the probability of a singly-occupied double dot. Such a 
pseudo spin representation is often used in the quantum 
information community. 4,5 With this change of variables, 
the dynamics of the double dot system can be mapped 
on the motion of a spin in an external magnetic field. 
This is in close analogy to the dynamics of a quantum 
dot connected to ferromagnetic leads. 28,29 

Due to the serial geometry the external tunneling af- 
fects only the z-direction of the pseudo spin and the left 
and right contacts couple with a different sign to I z . 
This is captured by the definitions £ll = (0,0,1) and 
riR = (0, 0,-1), which can be understood as pseudo-spin 
magnetizations of the leads. With this definitions the 
occupation probabilities obey the following master equa- 
tions: 



0=|Po = E %- gcT f r (E)P + \f-{B)Pi) 



(Bl) 



= J t P 2 = E Y ( f friE + U)Pl - f r( £ + U ^ 



Pl = l-P -P 2 



(B2) 
(B3) 



In equilibrium (/r = /l) the diagonal matrix ele- 
ments are given by the Boltzmann statistics Pq = 
l/Z,P l = 2g a cxp[~E/k B T}/Z, P 2 = g 2 a cxp[-(E + 
U)/k B T]/Z, Z = P + Pi + P 2 and the accumulation 
term as well as all components of the pseudo spin vanish. 

The dynamics of the single particle state is described 



Three different terms can be identified in the Bloch equa- 
tion. The term (dl / dt) acc describes the accumulation of 
pseudo spin in z-direction due to the serial external cou- 
pling. 

The relaxation-like term (dl/dt) rc \, limits the amount 
of pseudo spin. I relaxes isotropic by electrons leaving 
or entering the singly occupied double dot destroying all 
pseudo spin components. 

The third term looks like a rotation of the pseudo 
spin around a fictitious magnetic field B = (—A, 0, e re n), 
where e ren denotes the renormalizcd level separation 



= ( E L - Rc 



X 



(--) 



(E) 



g„xi + '-\E + U) 



-{g a - 1)X[ + -\E + U')\ ) — (L — ► R) , (B5) 

where Re denotes the real part. The Cauchy principal 
value integrals arc defined in Eq. (A2). 

This third term describes coherent oscillations inside 
the double dot which mix the accumulated spin in z- 
direction with the other components. The interdot tun- 
neling characterized by A leads to a precession of the 
isospin around the x-axes, while the energy separation 
between the dot levels results in a rotation around the z- 
axes. It is important to note that the renormalizcd level 
separation between the dots changes due to the external 
coupling and it is not given by the bare level separation 
e. 

In the following we outline the close analogy between 
the transport through a serial double dot described here 
and the dynamics in a spin valve described in Rcf. 28,29. 
Fig. 7 sketches a spin valve, realized by a single level 
quantum dot placed between anti-aligned ferromagnets. 
Relating the pseudo spin /, in the present work with the 
real spin S, in such a spin valve, one can perform the 
following mapping. The serial setup for the double dot 
system corresponds to the anti-aligned magnetization of 
the contacts in the spin valve. Furthermore the interdot 
tunneling translates to a transverse magnetic field in the 
single dot, while the level separation e corresponds to 
the magnetic field component along the magnetization of 
the contacts. Finally the renormalization of the energy 
levels discussed here was introduced in the spin valve as 
an exchange field leading to the Hanlc effect. 29 
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FIG. 7: Sketch of quantum dot spin valve. A single level quan- 
tum dot is connected to two ferromagnetic reservoirs with 
antiparallel magnetization. The spin precesses around an ex- 
ternal field with a component transverse (B±) and along (-By) 
the magnetization of the leads. Bu is modified by an exchange 
field arising due to the external coupling. This exchange field 
is manifest in the transport properties of the spin valve. 



According to equation Eq. (6) the stationary current 
is just given by the y— component of the isospin 



ieA 
~2h 



(P 



Pi) 



eA 

h y 



(B6) 



The system of master equations can be solved analyti- 
cally and the current as function of bias voltage and gate 
voltages has the following form: 



.4 



£ 2 + B 2 

°ren ' 



(B7) 



with the factors 

a Za . R 2 Z o , &? z b 
= lv~' = T + ^T 

Za = C^^rf rl +9a^rfr2) 
r 

N = 9a ^2 T r(fr2 + 9afr2){.frllfl + QafrlM) 
r 

+ T r(frl + 9afrl)(.f~ 2 fn + 9<jfr2ff 2 ) 
r 

{^9a{fL2 - /i?2)(^ r ''/'-l) 
HfLl-fm)CE r rfr2?j 

Zo ( (Er r r (/rl/r" 2 + Wrl/^ + tifrlfa)) 



Z A 



ZB - 4 \ T L T R 

+Jl2 (fm + Wm) + f R2 (fm + *9*hx) 

Here f denotes the opposite of r and we use the abbre- 
viations f rl = f r {E), f~ ± = 1 - / rl! f r2 = f r (E + U), 
f~ 2 = 1 — f r2l as well as the approximation f r (Ei,) w 

fr(E R ) P3 f r (E). 
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